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We give a numerical analysis of long- wavelength modes in the WKB approximation of cosmological 
scalar fields coupled to gravity via £4> 2 R. Massless fields are coupled conformally at £ = 1/6. 
Conformality can be preserved for fields of nonzero mass by shifting £. We discuss implications for 
density perturbations. 
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£ = g^d^d^ + ^ 2 R - V(</>) . (1) 



Q , I. INTRODUCTION 

O ; 

Scalar fields have played a major role in attempts to model the early Universe. In particular, nearly every incarnation 
O ■ of the inflation scenario has relied on scalars to generate vacuum energy and in turn exponential expansion and density 
fluctuations. Many of these models rely on slow roll potentials, i.e., potentials that are nearly flat where the scalar 
00 ' masses can be very small. 

Recently, horizon size and super horizon size density perturbations have been studied intensively, because of their 
, importance for understanding low i modes (see and references therein) in the WMAP data Qj. Long wavelength 
^ \ scalar field modes have interesting properties when the wavelength is on the order of the horizon size cHq 1 . One finds 
dispersion and diffraction effects that depend on the scalar mass and it's coupling to gravity. The generic lagrangian 
• for a scalar in a Friedmann-Robertson- Walker (FRW) Universe is 

o 
in 

If V((j>) contains no dimensionful parameters, then the scalar field is conformally coupled when £ = |. Conformal 
*-£h I invariance can be broken by including a mass term in V (</)). Here we assume the local (Minkowski limit) real tfi 4 theory 
is renormalizable. While this is not completely general, it is sufficient for our purposes. One could easily generalize 
our analysis to complex fields or fields in irreducible representations of some continuous symmetry group. 

One expects the scalars to be an integral part of any realistic model, so for instance, if the overarching theory 
is based on strings with a local or global SUSY preserved down to some scale, then the scalars will be components 
of some superfield $ contributing to the superpotential W($). This will put constraints on V(4>). In particular, 
flat directions could result (regions of moduli space where the scalar mass vanishes) and lead to massless or nearly 
massless modes, where for example SUSY could be broken by nonperturbative effects. We give these comments as a 
justification for the study of scalar zero modes and modes of very small positive mass or modes of very small imaginary 
mass (where the field could be rolling). However, there could be other reasons to study such modes. 

As the wavelength approaches the horizon size, the naive redshift formula no longer applies and one must refine the 
flat space analysis of the scalar field dispersion relation 0, Q ■ We will carry out a numerical analysis of the behavior 
of long wavelength scalar field modes and investigate the dependence of the redshift on the scalar field mass, and 
its coupling to gravity. We begin with a summary of scalar fields coupled to gravity in an FRW universe. We then 
review redshift and physical wavelength formulas, after which we are in position to begin our numerical analysis. We 
conclude with a discussion of the implications of our results. The interpretation and application of our results requires 
some comments. The Universe has expanded through a number of phases during its lifetime. We think we are now 
transitioning from a matter-dominated phase to a vacuum-dominated phase. Earlier there was radiation domination, 
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and before that, inflation, which occurred some time before Big Bang nucleosynthesis, but it is not clear how much 
before. Our task is to follow modes through these phases. 

It is unlikely that long wavelength modes can be measured directly so we are obliged to consider indirect measure- 
ments. These involve the long-wavelength background on which CMB or other short wavelength radiation propagates, 
and a proper analysis would consist of comparing observation with results predicted with and without dispersion. 

Our best guess for the relevant long wavelength modes that will lead to distortion of the CMB are those modes 
that were produced during inflation, then pushed outside the horizon during inflation, and have recently re-entered 
our horizon. These are the lowest I modes. They will have their distortion preserved due of the fact that they have 
spent the time from which they left the horizon until the present epoch frozen, and so unable to grow or dissipate. 
Being of the order of the present horizon size, they will also display the maximum distortion. We give generic results 
that can be applied to any model, but since results are inflationary model dependent, a full analysis would require 
specifying model parameters £, m, etc. 



II. SCALAR FIELDS COUPLED TO GRAVITY IN AN FRW UNIVERSE 

We first review the properties of a real free massive scalar fields coupled to gravity. The most general action (to 
first order in R) for this case is 

S = f d 4 xV=g{g^d^d„(f> - m 2 (f> 2 - ^ 2 R) , 



(2) 

where R is the scalar curvature and £ is a dimensionless coupling. We will work in FRW geometry, and use the 
convenient conformal parametrization of this family of spacetimes, 



9tiv = C(?7)diag(l, -/ly(x)) , 
where C(rj) = a 2 (t) is the conformal scale factor related to a conformal time via 



v(t) = 



1 cdt' 
o(f) 



The spacial part of the metric is 



ft-ij(x) = diag((l — nr ) ,r ,r sin 9) , 



(3) 



(4) 



(5) 



with k = 0,1,-1 corresponding to (flat), deSitter (positive) or Anti-deSitter (negative) curved spatial sections, 
respectively. 



The action J2J leads to the field equation 



□0 + m 2 (t> + £,4>R = . 



(6) 



where □ = -^i— gg^ v d v ). Because of the homogeneity and isotropy of the FRW metric, the solution to the field 
equation factorizes to 



Mv,x) = c- 1/2 ( v )h(v)yk(x), 



where 34 is an eigenfunction of the spatial laplacian 

1 



--di 



hh l 3d 3 y k (*) =-(|k| 2 ~ K )y fc (x), 



and k = |k|. In the massive case, the temporal part of has to satisfy 



fk + 



k 2 + m 2 C{7 1 )+£,--R{i 1 )C{ V ) 
6 



fk = 0. 



(7) 



(8) 



(9) 



3 



One can express the scalar curvature R in terms of scale factor C and curvature constant k in the form 



(10) 



Thanks to this relation, equation takes the form 



C 1 C 
C ~ 2 I C 



fk = 0. 



(11) 



The theory is conformally coupled if £ = jr and to 2 = 0. We will concentrate on two realistic cosmological regimes: 
vacuum (VDU) and matter (MDU) dominated epochs. In these cases (|llf) reduces to 



fk + 



(3 2 + m 2 C - 



fk = 0, 



where we have introduced the index v defined by 



^ 2 (£,p) = 7 -(6£-i) 



P(2p-1) 

(p-l)2 



(12) 



(13) 



(with p = 2/3 for MDU and is also formally 2/3 for VDU) and a conformal wave number j3, corresponding to a mode 
k: 



P 2 



4tt 2 

+ (6t-l)(n -l)H 2 



a 2 (t ). 



(14) 



Here Ao denotes the physical wavelength, corresponding to the wave number k, as measured today, f2o is the present 
ratio of matter-energy density to critical density and Hq is the present value of the Hubble parameter. Since all 
current observational evidence points toward a flat universe, we set n — (f2o = 1) so that (|14fl reduces to 



2ir 

(3 = k= —a(t ) ■ 



(15) 



A. Massless case 



In the massless case (jllf) reduces to 



,2 1 1 



ii 2 



fk = 0. 



(16) 



The solutions to this equation can be written in terms of Hankel functions He, (0rf) , which in polar form arc 

Hl 1 \prj)=A{pr,)e- iS ^ v \ (17) 

with A and S being real valued amplitude and phase functions. These are easily expressed in terms of ordinary Bessel 
functions J„ and Bessel functions of the second kind Y v . The phase is 



S{fir\, v) — arctan 



for real v(£,,p) and 



S(j3rj, v) = arctan 



cot(TTv)Jv{(3v) — csc{-Kv)J- v {(3rj) 

Jam 

[e-^JAPv) - J-uW)] 



%t[e- iv «J v {Pv)- J-uifiv)] ' 



(18) 



(19) 
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for imaginary v(£,p), while the amplitude is 



AWw) = y/J3(Pv)+Y*(J3r,), 

for real v and for imaginary v we find 

A(J3ri,u) = \csc(nri)\yJ{^[e- i ^J u (l3Ti) - J-uWv)]} 2 + S \e^MPv) - J-APv)}Y 
The angular frequency of FRW modes associated with a wave number k is given by 

dS 



drj 



where S is the corresponding phase given by either (|18Jl or 119(1 , depending on the choice of cosmology. 



(20) 



(21) 



(22) 



B. Massive case 



We want to write l|ll|) in the form fa + io\fk = 0. Therefore using a WKB analysis one finds the frequency (|22|l to 
second order fa + wj^ 2 /fc — 0, where uj^ is given .4] by 



4%) 



,(0) 



3£-± 



2w 



(o) 



C 
C 



c 
c 



(OK 



c-c- 



3m 2 C 2 
(o)> 



with 



,(°) 



a//? 2 + m 2 C{r 1 ) 



(23) 



(24) 



We have checked the validity of the WKB approximation (see appendix- lA"|l by comparing with the massless limit 
where exact solutions are available. The zeroth order contributions to the frequency and conformal scale factor and 
their derivatives for the cosmological cases of interest are summarized in Table-[IJ . 
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a° vac ~ 2.27 x 10 _18 s _1 
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-ca° ac (z + 1) 






MDU 


o .§ 

"mat*' 


aLat = [GttGnPo(w)\ w = o] 3 


a° m J -2.05 x lO" 18 ^ 1 


o 2 a W 
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TABLE I: Important factors for VDU and MDU. 

The scale factor and hence conformal time depends on the given epoch, so for power law expansions we have 

a(t) = [&kG n (1 + w) 2 Po (w)t 2 ] J ^ (25) 

as can be directly determined from the Friedmann equation, where w is the proportionality constant in the equation 
of state P(t) = wp(t) appropriate for a given background, and po{w) is a present value of critical density for a given 
epoch. Relevant choices of parameters for use in (|23|l and (|25l) are given in Table-[I]and Table-ITTl 



In the following we take the value of the present energy density to be po = Pcrit = 9.21 X 10 



-27 fcg 



4.42 X 10 17 s. 
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Epoch: (w,p) 



,(°) 



,(<>) 



a(t) 



c(v) 



Vacuum: 



(3 2 + m 2 



m 2 r] 3 



,(°) 



exp 



-po 



c 

To - 



c 

To - 



c 

To - 



Matter: [ 0, - 
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m 2 rf 



[6nG N p ] 3 t 3 



„0 6 
u mat 4 



o° 6 

, "mat 3 



81c 



TABLE II: Zeroth order result for frequency, scale factor and their derivatives, where a given epoch is characterized by: w, the 
proportionality constant in the equation of state P(t) = wp(t) appropriate for a given background, as well as the exponent of 
a power-law type cosmologies, i.e. a(t) ~ t v . 



C. Redshift formula 



The classical redshift formula in terms of frequency v is 

v a(t) 



(26) 



v a(t ) ' 

We want to find the correction factor to this naive redshift formula where the correction is the result of the nontrivial 
modifications to the dispersion relations for long- wavelength modes i.e., wavelengths of the order of the horizon size. 

• n 

To do this we have to take into account of the conformal angular frequency correction 3 so that we find 

= a(t) cj k (t a ) 
v a(t ) Lu k {t) 

Here for uok we use cu k ^ = {$ in the massless case and u k 2 \ given by equation (|23|l . in the massive case. In the 
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t k 
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Xphy = a(t)\ 


frequency 


_ dS 
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LU 



TABLE III: Dictionary of conformal and physical variables. 



following sections we are going to present the results for two different cosmologies. The relations between physical 
and conformal variables are summarized in Tables [I] and IIIII The advantage of an analysis via the WKB method 
is that it is simple and straightforward, and it usually gives the correct trends when the corrections are moderate 
(~ 5% to ~ 20%) (as we have shown in AppcndixjXJ . These observations can be verified by comparing with exact 
results where they exist. In the cosmological regimes of relevance (vacuum and matter domination), we find the 
dependence of the dispersion relations on the value of the mass of the scalar field. We formulate our discussion in 
terms of conformal invariancc, i.e., in terms of the value of the mass and wavelength where the conformal behavior is 
approximately preserved. 

In all cases we find m ~ 10 _33 ey (inverse Hubble size) as the mass where nonconformal behavior starts to set in. 
These masses should be nearly equal in the different regimes, since differences are caused by numerical factors of order 
one. 



6 



To proceed further, we first have to express all the parameters present in eqs. (|24fl and i|23|) in terms of redshifts 

a(to) 



a(t) 

and the parameter bo 



(28) 



feo^^r, (29) 
cH 

that can be interpreted as the fractional size of the physical wavelength in the units of the present Hubble radius. 
We will work in units where a(t) and hence C(rj) are dimensionless. The wave number [3 does not depend on the 
epoch, as can be seen from equation (|15(l . We have now collected all the necessary epoch specific input needed for 
our numerical analysis that will be carried out in the next section. 



III. CORRECTION TO REDSHIFTS 



In this section we are going to present corrections due to the redshift formula Q27JI originating from the coupling of 
the scalar fields to gravity. We consider massless as well as massive cases in two different cosmologies, i.e., vacuum and 
matter dominated universes. In order to see the full spectrum of possible behavior of the dispersion relations, one has 
to discuss both real and imaginary masses. We plot the ratio of initial to final frequencies in each epoch. Sequentially 
through matter and vacuum domination, we set z™ at — 1100, next z™ at = z^ ac = 10 and finally zj ac — 0. In the case 
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FIG. 1: Vacuum Domination m 2 > WKB: Dashed curves £ < 1/6, Thick curves £ > 1/6. In this and the following figures we 
use £ = ±3/4, ±1/2, ±1/4, for fs. 



of vacuum domination, for real masses (Fig-^l, the deviation from the classical redshift formula is very small until 
length scales of the size ~ 0.6c Hq . For imaginary mass, when 9(m) < 10~ 32 eV, the ratio of frequencies exhibits 
similar behavior (Fig-|2J). In both cases, as the magnitude of the mass of a scalar field gets larger than roughly 10~ 32 
eV, it dominates the effects of £ if £ < 1. In the matter dominated case, both real (Fig-|3J), and imaginary (Fig-QJ 
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FIG. 2: Vacuum Domination m 2 < WKB: Dashed curves are for £ < 1/6, and thick curves are for £ > 1/6. 



masses of order ~ 1 10 31 1 eV^ dominate effects when £ < 1. The corrections become substantial at smaller length scales 
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FIG. 3: Matter Domination m 2 > WKB: Dashed curves are for £ < 1/6, and thick curves are for £ > 1/6. On this and the 
following figures vertical lines show the corresponding asymptotes where the WKB approximation fails. 
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FIG. 4: Matter Domination m 2 < WKB: Dashed curves are for £ < 1/6, and thick curves are for £ > 1/6. 



IV. WMAP FIT AND DISCUSSION 



Now let us discuss the processing of the density perturbation spectrum. As is well known, once a perturbation 
comes within the horizon, it begins to oscillate. This phenomenon is reflected in the observed large £ WMAP CMB 
spectrum, where the first maximum (first acoustic peak) is at I ~ 200, and higher order peaks are at larger I value. 
The low £ values have not been inside the horizon long enough for much processing to have taken place. The region 
from roughly £ — 20 to £ — 200 is just now beginning to undergo its first plasma oscillation, while for £ < 20 very 
little has happened yet. But this is just the region of interest for the dispersive effects we have been discussing, and 
we are lucky that in this region (£ < 20) we have a pristine, unprocessed spectrum. (Recall for large £, we have short 
wavelengths, and so virtually no dispersion.) Hence we can confine ourselves to the analysis of perturbations with 
wavelengths of the order of or somewhat less than the horizon size where we do not need to worry about plasma 
oscillations. (There may be issues of re-ionization to consider, but at the level we are working we choose to ignore 
such effects.) 

Perturbations that are of the order of the horizon size today have undergone an evolution from the time of their 
production. A typical scenario would be: Perturbations are produced during a vacuum-dominated epoch of inflation, 
and are pushed outside the horizon. Inflation then ends, and the universe becomes radiation dominated. Some 
perturbations come back inside the horizon, are processed via plasma oscillations, etc., until about z = 1100, when 
the universe becomes matter dominated. More perturbations re-enter the horizon and are processed until around 
z = 10 when the universe again becomes vacuum dominated and perturbations again start to be pushed outside the 
horizon. Ultimately, we would like to follow the entire evolution of a perturbation from its production until today for 
modes that are currently near horizon size (an £ < 20 mode). However, this would require a detailed model of the early 
universe. A less ambitious approach is to follow some (better understood) fraction of the evolution to demonstrate 
that dispersion can play a role in understanding the observational data, and leave it to future work, when a more 
detailed understanding of the early Universe including details of early Universe phase transitions are known, to follow 
the complete evolution of the modes. To this end we have shown in Figure JSJ (where we convert wavenumber k 
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to multi-pole moment £ using k = £/(cHq 1 )) the results of evolving modes from 2 z = 1100 until today (i.e. from 
z = 1100 to z = 10 with matter domination, and from z = 10 until today (z = 0) with vacuum domination), and 
have fit the results to the low £ WMAP data. Since an £ = 20 mode undergoes very little processing or dispersive 
evolution, we have normalized our amplitude to this region of the spectrum. Once this is done, we have a single free 
parameter £ for massless scalar fields (p. We then do a one-parameter fit, as shown in Fig-0 and find £ ~ I — 0.0002 
see Fig-|51 This is very close to the conformal value £ = | and can be shifted there, but only at the expense of of 
introducing a small negative mass squared for the scalar field as can be seen from the field equation (fill) evaluated 
for the VDU 



fk{v) 



2 , (m-2) + (mc/a° vac f 



It is clear that one can introduce the effective coupling £ such that 

12f = 12f - 



fk = . (30) 



(31) 



and by setting £ to the conformal value 1/6 we can find the value of mass m corresponding to a field <p with an 
effective conformal coupling to gravity. Our conclusion is that the evolution of large £ scalar modes can be used to 
constrain the coupling of scalar fields, in particular the inflaton, to gravity. If we set m — we find that minimally 
coupled field (£ = 0) is easily excluded. Our fit is merely an example of how constraints on £ can arise. Specific models 
will give specific results. It is interesting that, with a few assumptions, a value of £ can in principle be extracted 
from a study of the cosmic microwave background. The coupling of scalar fields to gravity have other ramifications 
that would need to be considered in any realistic model. Another reason for being cautious about drawing sweeping 
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FIG. 5: Fit of the low I part of WMAP spectrum to scalar fields with dispersion. The thick line is for the best fit value, 
£ = 0.166435. We have added three other curves for comparison. The upper thin line is for £ = 0.169492, the flat line is for 
£ = 1/6, and the lower thin line is for £ = 0.161290. 

inferences from Fig-|S]is that the action given in with m 2 = leads to a spectral index in disagreement with the 
data. Perhaps what one should conclude is that we need a theory with a sufficiently complicated potential V((j>) 
that density perturbations can be laid down when the <f> mass is sufficiently large (see |16J), but where the late time 
effective theory is nearly conformally invariant. 
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FIG. 6: Deviation of £ from the conformal value of coupling £ = i as seen on a plot of x 2 vs - (§ — £,)■ 




FIG. 7: Percent error in VDU (a) and MDU (b) for both negative and positive values of £. 



APPENDIX A: VALIDITY OF WKB APPROXIMATION 



The WKB approximation works best for couplings close to the conformal case £ = 1/6. For various choices of 
parameters, we show the percent errors of the WKB approximation relative to the exact vacuum dominated era 
results (Fig- 7(a) and Fig- 8(a) I, and to the matter dominated universe in (Fig frnj}| and Fig- p^b)} . In a matter 
domination universe represented in Fig- 8(b) the WKB approximation is better (i.e. up to larger scales) for positive 
couplings. However, the broad range of applicability of the WKB method allows us to use the WKB approximation 
to draw conclusions about trends in the data. 
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FIG. 8: Comparison of exact solutions to the WKB approximation in VDU (a) and MDU (b) case. Thin (thick) lines represent 
exact (WKB) solutions with £ > 1/6. Small (large) dashing represents exact (WKB) with £ < 1/6. 
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